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EXTRAGALACTIC SETI: THE TULLY-FISHER RELATION AS A PROBE OF DYSONIAN 
ASTROENGINEERING IN DISK GALAXIES 

Erik Zackrisson 1,2 *, Per Calissendorff 2 , Saghar Asadi 2 & Anders Nyholm 2 

ABSTRACT 

If advanced extraterrestrial civilizations choose to construct vast numbers of Dyson spheres to har¬ 
vest radiation energy, this could affect the characteristics of their host galaxies. Potential signatures 
of such astroengineering projects include reduced optical luminosity, booste d infrared luminosity and 
morphological anomalies. Here, we apply a technique pioneered by [Anni 1 (119991) to search for Kar- 
dashev type III civilizations in disk galaxies, based on the predicted offset of these galaxies from the 
optical Tully-Fisher relation. By analyzing a sample of 1359 disk galaxies, we are able to set a con¬ 
servative upper limit of < 3% on the fraction of local disks subject to Dysonian astroengineering on 
galaxy-wide scales. However, the available data suggests that a small subset of disk galaxies actually 
may be underluminous with respect to the Tully-Fisher relation in the way expected for Kardashev 
type III objects. Based on the optical morphologies and infrared-to-optical luminosity ratios of such 
galaxies in our sample, we conclude that none of them stand out as strong Kardashev type III can¬ 
didates and that their inferred properties likely have mundane explanations. This allows us to set a 
tentative upper limit at < 0.3% on the fraction of Karashev type III disk galaxies in the local Universe. 
Subject headings: Extraterrestrial intelligence - galaxies: spiral - galaxies: stellar content infrared: 
galaxies 


1. INTRODUCTION 

For more than 50 years, astronomers have searched the 
skies for artificial sig nals and oth er sign s of extraterres¬ 
trial civilizations (e.g. lTartei!200lLl2007lh yet no evidence 
of intelligent life beyond the Earth has so far emerged. 
One viable explanation, among the many pro posed f or 
this so-called Fermi paradox (e.g. lHrTnllT983l: fWebhl[200l 
lGirkovi(jl2009 ) is that technologically advanced civiliza¬ 
tions are very rare. Most SETI efforts have so far been 
limited to the confines of the Milky Way, and with > 10 11 
large galaxies in the observable Universe ([Beckwith et al.1 
2006:), it has been argued that the prospects of detection 
may be better if the search radius is extended to ex¬ 
tragalactic scales (e.g. iCirkovic fe Bradburvll2006t iWilevl 
1201 If) . The transmission of signals across extragalactic 
distances would be far more ch alleng ing than across a sin¬ 
gle galaxy (see e.g. lMaccondl2012l for a proposed tech¬ 
nique of this type), but extremel y advanced sup ercivi- 
lizations (Kardashev type II-III; iKardashevlll964 ) could 
in principle be detectable through indirect signatures of 
galactic-scale astroengineering projects even if they have 
no interest in intentionally transmitting signals in our di- 


rection (Annis 1999 

; Carrigan 201(1 ] 

noue & Yokoo 2011); 

Wright et al. 2014a 

.b; Griffith et al. 

20lJ). 


One potential way to harvest the radiation energy of 
a sta r is the construction of a Dyson sphere dDvsonl 
Il960f) . which would either completely or partially en¬ 
shroud the star and make a significant fraction of ra¬ 
diation energy available for supercomputing or ot her 
project with immense power requirements fe.g. lSandberd 
1 1 999f iBradburvl 120011 ). This would render the star un¬ 
usually faint or completely dark at ultraviolet/optical 
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wavelengths, but could make it bright at infrared wave¬ 
lengths due to the need to get rid of excess heat. Sev¬ 
eral unsuccessful searches for individual Dyson spheres 
within the Milky Way hav e already been carried out 
(I Slvshl 1 198 5t lTimofeevfl20dol : l.lugaku fc Nishimural 12004 
lCarrigadl2009f) . 

In a pioneering study, lAnnisI (|1999D proposed the use 
of galaxy scaling relations to single out potential Kar¬ 
dashev type III candidates, i.e. galaxies in which an 
extremely advanced civilization has spread throughout 
its host galaxy and converted vast numbers of stars into 
Dyson spheres. As such galaxy-scale astroengineering 
projects would likely not affect the gravitational poten¬ 
tial of the galaxy, but substantially decrease the total 
optical luminosity, host galaxies of Kardashev type III 
civilizations may appear anomalously underluminous for 
their type. In a sa mple of 57 disk galaxies and 106 ellipti¬ 
cals, [Ambl (U999H de t ected no significant outliers. Here, 
we revisit the lAnnhj (|1999H method and apply it to a 
Tully-Fisher sample of 1359 disk galaxies. 

The search method, the detected Tully-Fisher outliers 
and the resulting upper limits on the prevalence of star- 
fed Kardashev type III civilizations in disk galaxies are 
discussed in Sect. [2] In Sect. [U the nature of the outliers 
is probed further using morphological information and 
infrared data. Other potential signatures of supercivi¬ 
lizations that make use of Dyson spheres on large scales 
are discussed in Sect. [4] Sect. [5] summarizes our findings. 

2. KARDASHEV TYPE III CANDIDATES IN THE 
TULLY-FISHER DIAGRAM 

The Tully-Fisher relation (jTullv fc Fished fl977 l is an 
empirical relation between the velocity width (amplitude 
of the rotation curve) of a spiral galaxy and its absolute 
magnitude. Because of the low scatter in this relation, 
it is often used as a distance indicator: By combining 
an apparent magnitude derived from photometry with a 
spectroscopic measurement of the velocity width (usu- 
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ally based on the HI 21 cm line or the Ha line at 6563 
A), it becomes possible to infer the distance by adopt¬ 
ing the absolute magnitude predicted by the relation. 
The Tully-Fisher relation has been shown to hold for ro- 
tationally supported galaxies over an imp ressi on range 
in surface brightness and luminosity (e.g. iZwaan et al.1 

ll995tlMcGaugh et, al.lll998l : iChung et al.ll2002H . although 

there are indications of increased scatter and a system¬ 
atic change in slope for the faintest and s mallest galaxies 
([Matthews et al.lll998HGeha et al.ll2006D . 

As pointed out by|Anni si (I1999T) . galaxy-wide coloniza¬ 
tion projects in which large numbers of stars are con¬ 
verted into Dyson spheres would affect the absolute mag¬ 
nitude in short-wavelength passbands, but are not likely 
to affect the velocity width, which serves as a proxy 
for the depth of the gravitational well. The host galax¬ 
ies of star-fed Kardashev type III (hereafter Kill) civi¬ 
lizations should therefore appear significantly underlumi- 
nous compared to the optical Tully-Fisher (hereafter TF) 
relation. In search of Kill host galaxy candidates. lAnnia 
( 1999 1 adopted the criterion that such objects should 
have their optical luminosity dimmed by at least a fac¬ 
tor of four (i.e. 1.5 magnitudes) compared to the TF 
relation, and an alyzed a sample of 57 disk galaxies from 
iPierce fc Tullvl ( 1988 1 without finding any TF outliers 
according to this criterion. Here, we apply the same tech¬ 
nique to TF data (luminosity in the /-band and velocity 
width Wtf based on eit her the HI 21 cm line or Ha) 
for 1359 objects from the lSpringob et ahl (|2007l ) SFIH—|- 
catalog. 

The criterion that Kill galaxy candidates should be 
located > 1.5 mag below the TF is simply chosen for 
conv enience - the TF r elation has an intrinsic scatter 
(e.g. iMasters et al.ll20d6h that, when coupled to observa¬ 
tional errors would lead to too many spurious detections 
if the li mit had been placed clo ser to the TF. The bary- 
onic TF dMcGaugh et al.ll2000 j exhibits smaller intrinsic 
scatter and could potentially be used to search for less ex¬ 
treme Kill objects, but this requires additional data for 
each object (an estimate of the molecular gas fraction), 
and current baryonic TF samples are consequently much 
smaller. Objects located significantly more than 1.5 mag 
below the TF become increasingly likely to be excluded 
from existing TF samples due to selection effects (which 
we attempt to quantify in Sect. ESI). For these reasons, 
we have chosen to stick to the >1.5 mag criterion used 
by Anni si (1 19991). Using the PARSEC vl.2S isochrone s 
(|Bressan et al.ll20liilChen et al.ll20l4 iTang et aLll2014T) . 

as implemented in the CME0 spectral synthesis code, we 
estimate that 40-50% of the stellar contribution to the 
/-band flux of galaxies that have undergone active star 
formation for ~ 10 Gyr comes from main-sequence stars, 
whereas the remaining fraction comes from more evolved 
stars (red giants, horizontal branch stars, asymptotic gi¬ 
ant stars). Hence, a very wide range of stars would need 
to be enshrouded in close-to-complete Dyson spheres in 
order for Kill civilizations to be detectable with the TF 
method. 

Below, we present a brief overview of how the TF data 
in the SFI-I—I- catalog was originally derived. Further 
pruning and modifications of this catalog for the pur- 
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poses of our Kill search are described in Sect. 12.11 
In the iSpringob et all (I2007T) data set, Wtf is either 
derived from spectroscopy of the HI 21 cm line or from 
optical Ha rotation curves. In the HI case (the most 
common case for the TF outliers we find), Springob et al 
derive the line width Whi using: 


Wtf 


Whi. 


obs 


- A iri 


1 + 2 



1 

sini 


(i) 


where Whi, obs represents the observed line width of the 
HI 21 cm line, z the redshift (based on the wavelength of 
either the HI 21 cm line or the wavelength of Ha), i the 
inclination, Aj ns the instrumental correction to the line 
width and At the corresponding turbulence correction. 
In the case of A t , a value of 6.5 km s _1 was adopted. 
For a detailed description of the i n strum ental line width 
correction A t , see ISpringob et all (120051) . 

In the case when optical rotation curves are used in 
SFI++, corrections are applied for known systematic dif¬ 
ferences between Ha and HI velocity widths due to the 
limited extent of the Ha data. The inclinations used 
to correct Wtf for orientation effects in the catalog is 
derived using 


i = cos 


( ( b / a ) 2 ~ 9 o 

V 1 - <7o 


0.5 


( 2 ) 


Here, a represents the senhmajor axis, b the senhnhnor 
axis, qo the intrinsic axial ratio b/a of the disk at edge-on 
orientation, assumed to be qo = 0.13 for disks of morpho¬ 
logical type Sbc and later, and 0.20 for disks of earlier 
morphological types. 

The /-band absolute magnitudes of SFI++ have been 
corrected for internal dust attenuation by an amount 
AM att given by 

AMatt = - 7 log 10 (a/&), (3) 

where 7 is taken to be a function of Mj (with high- 
luminosity disks assumed to disp lay higher dust atten¬ 
uation; see ISpringob et ail 12007 for details). The data 
have also be en corrected for Ga lactic extinction using 
COBE data (|Schlegel et al .1 IT998li . 

Pinning down the factors that dominate the error bud¬ 
get on the position of an object in the TF diagram is 
highly nontrivial, given the many measu r ement s and cor¬ 
rections involved, but IGiovanelli et al.l (119971) attempt 
such an analysis for a subset of the SFI++ data. They 
find that the errors and corrections on the velocity width 
tend to dominate, except in the case of high-width, high- 
luminosity, highly inclined galaxies, for which errors on 
the extinction correction (eq. [3]) can become the main 
source of uncertainty. The median correction for dust 
in our sample is ~ 0.3 mag, but in extreme cases, the 
correction does reach ss 1.5 mag. Errors on photometry, 
redshift and ellipticity are typically negligble by compar¬ 
ison. Of course, this holds for normal galaxies, absent 
Kill activity. The sensitivity of the dust correction to 
the optical axial ratio (b/a) could in principle jeopardize 
the assumption that Dysonian astroengineering should 
affect the optical stellar luminosity (as measured by A//), 
while leaving the overall dynamics and the gas phase (as 
measured by Wtf) untouched. If Kill activity results in 
morphological curiosities which affect the apparent (b/a) 
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ratio, then this could affect both Mj and Wtf, through 
the inclination-related corrections going into both eq. m 
and eq. d3j) . 

Based on the numerical simulations of Kill coloniza¬ 
tion presented in Sect. 13.11 we argue that this effect can 
go both ways, and act to either accentuate or mask Kill 
activity in the TF diagram. An evenly colonized Kill 
disk with uniformly diminished surface brightness may 
appear rounder than it actually is (due to a spuriously 
high b/a ), since only the inner, bulge-dominated parts 
may remain visible above the surface brightness thresh¬ 
old in optical images. Such objects would be pushed 
towards the lower right of the TF diagram (deeper into 
the Kill candidate region), since Wtf would be biased 
high and the dust correction on Mj would be underesti¬ 
mated. An already highly inclined galaxy in which Kill 
colonization has focused on the bulge-dominated inner 
regions may on the other hand appear flatter than it ac¬ 
tually is (due to a spuriously low b/a), thereby pushing 
it upwards (towards the TF line and away from the Kill 
region). This sensitivity to how Kill colonization has 
progressed within a galaxy is an inherent limitation of 
the TF approach to Kill searches, as long as the axial 
ratio is based on the optical appearance, and should be 
kept in mind as a possible caveat when considering the 
TF-based limits on astroengineered disks we derive in 
Sect. 12.31 


2.1. Distance estimates 

In the original SFI+-f- catalog, distances (in velocity 
units) to galaxies not belonging to any group or cluster 
are estimated using the relation 


T’gal — CZgal ^gab 


(4) 


where cz ga i is the velocity in the CMB frame and u ga i the 
peculiar velocity expressed as 

Ugal = C2T gal (l - 10 O 2dm ). (5) 

Here, dm represents the difference between the 
attenuation-corrected absolute magnitude (as presented 
in the SFI+-f catalog) and the predicted absolute mag¬ 
nitude one would expect if the object had no offset from 
the TF relation. 

However, since we are specifically looking for TF out¬ 
liers, distances based on the assumption that a single 
object lies on the TF are of little use. For objects be¬ 
longing to galaxy groups or clusters containing several 
disk galaxies with individual TF distances, the average 
distance to the group or cluster itself may still be used 
as an estimate of the distance to the member galaxies. 

In our analysis, we therefore discard all field galaxies 
from the sample and estimate new absolute magnitudes 
Mi for the remaining objects based on the SFI++ dis¬ 
tances assigned to the groups and clusters to which they 
are deemed to belon g. _ 

In iSpringob et al.l (120071 ), such group and cluster dis¬ 
tances are estimated similar to eq. 0, but with cz gIp 
calculated as the mean redshift of all group members 
and peculiar velocities estimated as 


E*= 


i=l,JVsFI ++ Vl 


A 


i^grp 


y i/ e . ’ 

Z^i=l,AT SPI++ 1 / t l 
where Vi is calculated similar to eq. ©, with the slight 


difference of using cz glp instead of cz ga \ for individual 
gala xies. Here e, represents the estimated error on Vi. 
ISpringob et all (120071) also present a version of u grp cor¬ 
rected for Malmqvist bias, r grPj M (distance in velocity 
units), and this is the quantity we use to recompute the 
absolute magnitudes: 


Mi = m/c — 5 log 10 (r grPjM /-ffo) + 5. (7) 

Here, m/ c is the /-band apparent magnitude cor¬ 
rected for both internal and Galactic extinction and 
Hq is the Hubble constant. To be consistent with 
the luminosi ty-dependent ex ti nction correction scheme 
derived bv IGiovanelli et, all (119951 ) and adopted by 
ISpringob et al.l (120071 ). we have adopted H 0 = 100 km 
s _1 Mpc -1 when computing the extinction corrections, 
but have then rescaled all absolute magnitudes to Hq = 
70 km s -1 Mpc^ 1 , which is more in line with modern 
measurements of the Hubble parameter. One should 
note that the extincti on correction scheme derived by 
IGiovanelli et all (Il995l) was based mainly on Sbc and Sc 
galaxies, whereas SFI++ contain s a gre ater diversity of 
spiral subclasses. ISpringob et al. ( 200 7) argue that, be¬ 
cause of this, extinction in earlier types may possibly be 
underestimated. However, this turns out to be of lit¬ 
tle consequence for the present study, since very few Sa 
galaxies end up as Kill candidates (see Sect. [3] for more 
discussion on Sa subtypes). 

The underlying assumption when giving preference to 
group distances over the distances to individual objects 
given in SFI+- 1- is that star-fed Kill civilizations are re l¬ 
ative ly ra re (as already indicat ed by the studies of lAnnisI 
119991 and iWright et al.l l2014bl) . so that even though a 
single group/cluster member may be a TF outlier, the 
TF distance to the group/cluster as a whole need not be 
significantly biased. In the case of a colonization wave 
spreading across intergalactic distances from one galaxy 
to the next, this assumption may of course be jeopar¬ 
dized, as further discussed in Sect. [4] A second caveat 
is that, since groups/clusters can have non-negligible ex¬ 
tent along the line of sight, the distance to a member 
galaxy may be somewhat different from the distance to 
the group/cluster centre. 

After rejecting field galaxies from SFI+- )-, we further 
prune the sample by removing galaxies belonging to 
groups and clusters with less than 3 SFI++ members, 
since too large distance errors would otherwise be in- 
trouced (see discussion in Sect. 12.21) . Out of 4861 galax¬ 
ies in the full SFI+-(- catalog, our remaining subsample 
consists of 1359 objects. Our scheme for reassessing ab¬ 
solute magnitudes for group and cluster members then 
generates a sample of 11 galaxies seemingly located >1.5 
mag below the TF relation, which th ereby qualify as Kill 
candidates according to the I Annisl ( 1999 ) criterion. The 
properties of these galaxies are summarized in Table [1] 

As it turns out, our 11 Kill candidates all belong to 
the so-called nontemplate groups of the SFI++ cata¬ 
log, i.e. groups o f galaxies with rela ti vely few other 
SFI+ + members ({Masters et all 120061 : ISpringob et al.1 
120071). As the exa c t spa tial sizes of the groups defined 
bv ISpringob et all (120071 ) are not known, we introduce 
an uncertainty on the absolute magnitudes by assigning 
each group an estimated radius of 1.5 Mpc. Assuming 
this to be the outer edge that confines the group, we 
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calculate the absolute magnitude that our objects would 
have at this distance and adopt the difference between 
the central and edge magnitude as the uncertainty. As 
a result, the uncertainty in absolute magnitude conse¬ 
quently grows with decreasing group distance (reaching 
a maximum of ±0.25 mag for the most nearby object in 
Table [T]). In order to include uncertainties from extinc¬ 
tion and inclination we also add the error of th e absolute 
magnitude tabulated by ISpringob et al.l (l2007h to the to¬ 
tal error budget. 

2.2. Quality of the Kill candidates 

One of the risks involved in deriving absolute magni¬ 
tudes for galaxies based on the TF distances to their 
parent groups, is that this may result in spurious Kill 
candidates in situations where background galaxies have 
been incorrectly assigned to foreground groups. This 
leads to underestimated galaxy distances and therefore 
absolute magnitudes that are much too faint. A poten¬ 
tial indication of this is when the redshifts of faint out¬ 
liers appear to be significantly higher than those of the 
groups hosting them. Peculiar motions may account for 
some of these offsets, but for four cases in Table Q] (ob¬ 
jects 1, 2, 9 and 10), the differences between the observed 
galaxy redshifts cz D b s and the group distances r grPi M cor¬ 
respond to velocity offsets of more than 3000 km s -1 , 
which seems extremely high. However, scrutiny of the 
SFI++ catalog entries reveals that the observed redshifts 
of these four galaxies agree well with the average group 
redshifts, and that it is the group TF distances r grPi M 
themselves that are dubious. These four outliers belong 
to two small groups with only three TF members each. 
In both cases, the third member happens to be a low- 
mass galaxy (log 10 Wtf ~ 2) with very high luminosity 
(« 3 magnitudes above the TF relation). Hence, the in¬ 
ferred TF distances to these groups seem to be biased by 
small number statistics and TF outliers, albeit outliers 
in the opposite sense of what we are looking for. 

Since the adopted distances to the groups hosting our 
Kill candidates are sensitive to the number of other TF 
galaxies in these structures, groups with many TF mem¬ 
bers are in general expected to have more reliable dis¬ 
tance estimates than those with few members. The 11 
outliers in Table Q] are part of 7 groups, where the small¬ 
est ones are the groups with SFI++ identifiers 40702, 
30590 and 31043 (3 members each in the SFI++ cata¬ 
log) and group 40032, (4 members). The largest groups 
are group 30654 with 29 members, group 1736 with 10 
and group 3407 with a total of 8 members. Here, we 
have assigned a default class ’B’ to candidates with only 
3 TF members, as these are more prone to errors, and 
class ’A’ to groups with > 4 members. Because of the 
clear-cut problems with the groups hosting objects 1, 2, 
9 and 10, these are further degraded from class ’B’ to ’C’, 
and are deemed to be the least likely Kill candidates in 
Table [L] When comparing the offsets between observed 
galaxy redshifts and the mean redshifts of their parent 
groups, objects 5 and 8 also stand out, with velocity off¬ 
sets of « 500 km s -1 and ss 700 km s -1 respectively. 
Since these are members of the two largest groups, with 
SFI++ velocity dispersions similar to or larger than these 
offsets, this is by itself not a source of concern. However, 
in the case of object 8, which would belong to class A 
according to our criteria, we suspect that there is some- 



Fig. 1. — Tul l y-Fish er diagram containing the 1359 objects from 
ISpringob et al.l (120071) discussed in Sect. [2] with error bars sup¬ 
pressed for all objects except the Kill candidates to avoid clutter¬ 
ing. The solid line indicates the TF relation and the dashed line 
the li mit below wh ich outliers qualify as Kill candidates according 
to the lAnmsl d 19991) criterion. Numbered dots mark the 11 outliers 
discussed in the text, with blue indicating the Kill candidiates of 
class ’A’, green indicating class ’B’ and red indicating class ’C’. 
Some of the errorbars have been slightly exaggerated for plotting 
purposes. 

thing; wrong; w i th the data or the catalog entry, since 
ISpringob et al.l (I2007T ) lists an apparent magnitude cor¬ 
rected for inclination and extinction that is fainter than 
the uncorrected observed magnitude, which is unphysi¬ 
cal. This objects has therefore also been transfered to 
class ’B’. 

Outliers of class A, B and C are marked in Fig. [T] as 
blue, green and red symbols respectively. As seen, the 
outliers with the least reliable distance estimates (class 
C) are the ones with the most extreme offsets from the 
TF relation, whereas objects belonging to class A and B 
straddle the Kill delimiter line. 

It should be noted that the present analysis is in¬ 
sensitive to extreme Kill objects, for which > 90% of 
the optical light has been absorbed by Dyson spheres 
(i.e. galaxies located > 2.5 mag below the TF), since 
such disks would not make it into the SFI++ sample in 
the first place. A number of such “dark galaxy” can¬ 
didates, with very faint or undetected optical counter¬ 
parts, have_already_been uncovered in blind HI sureys 
(e.g. lHavnes et al1l201lH . Most of these ha ve ve ry low 
masses and velocity w idths (e.g. I Cannon et al.l 1 201 51 
IJanowiecki et al.l 120151) . and are therefore in the region 
where t he TF a lrea dy is known to exhibit considerable 
scatter (|Geha et al.ll2006l l. Whether there are massive 
dark galaxies remains an o pen question - while such de- 
tections have been claimed (|Minchin et al.lf2005l 12007T) . it 
has also been argued that these could simply be tidal HI 
structure s with extended HI line profiles without ordered 
rotation (iDuc fe Bournaudll2008D . 

2.3. The prevalence of Kardashev type III civilizations 

It is important to realize that a substantial fraction 
of galaxies in the SFI++ compilation would likely have 
been excluded from the sample, had they been converted 
into Kill objects, as this would have rendered them too 
faint for detection or otherwise altered their properties 
in ways that would have prevented them from entering 
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the catalog. In fact, the prospects of finding Kill can¬ 
didates quickly deteriorate as the condition of TF offset 
A (Mj) is sharpened. The method becomes far less effi¬ 
cient at A (Mi) = 2.5 mag (90% of light absorbed) than 
for A (Mi) = 1.5 mag (75% of light absorbed; our de¬ 
fault limit), and fails to produce any useful results at 
A (Mi) = 5.0 mag (99% of light absorbed). 

To convert our harvest of seemingly underluminous 
disk galaxies into a constraint on the prevalence of star- 
fed Kill civilizations in disk galaxies, we need to quan¬ 
tify the fraction of objects that would not have been de¬ 
tectable at a given A (Mi). However, this become highly 
non-trivial becau se o f the inhomo gene ous selection cri¬ 
teria used in the iSpringob et all (120071 1 compilation. A 
casual upper limit can nonetheless be derived from the 
effective surface brightness limit of the SFI++ sample. 
By considering the size-flux distribution of objects in the 
catalog, we find that Kill-colonized galaxies that end up 
with a mean /-band surface brigthness within the outer¬ 
most detected /-band isophote /Jj > 23.5 mag arcsec -2 
are likely to be too faint for detection. To estimate the 
fraction of Kill candidates that would have been ex¬ 
cluded for this reason, we assume that Kill coloniza¬ 
tion causes a uniform dimming of the disk, and hence 
a rescaling of the entire surface brightness profile while 
leaving its overall slope unaltered (see Sect. I3.1l for a dis¬ 
cussion on this). In this case, Kill colonization makes /// 
fainter, but also makes the galaxy appear smaller, since 
the faintest detectable isophote is reached at a smaller 
distance from the centre. The latter effect is important, 
since some of the subsamples of the SFI++ compilation 
are based on optical size constraints, which may reject 
small Kill galaxies even if they are sufficiently bright to 
be observed. 

To simulate these selection effects, we fit exponential 
disk profiles to all of the 1359 objects in the SFI+-(- cat¬ 
alog for which the relevant photometric data is avail¬ 
able (about 80% of the objects). We then lower the sur¬ 
face brightness uniformly by a constant factor (A(M/)), 
derive the radius of the outermost detectable isophote 
v max,/ and the mean surface brightness within this 
radius. Based on the apparent surface brightness limits 
of the sample, and strict size limits for some of the sub¬ 
samples in the SFI++ compilation, we consider object 
likely to be rejected from the sample if /!/ > 23.5 mag 
arcsec -2 or r maXj / < 30 arcsec. Through this procedure, 
we find that « 30%, or « 400 objects would likely have 
remained detectable after 1.5 mag of dimming (75% of 
the light lost). Given 11 candidates out of a sample of 
400, one may then conservatively conclude that < 3% 
of local disk galaxies are housing star-fed Kill civiliza¬ 
tions. These estimates also reveals the limitations of the 
method - only 16%, 9% and 3% of the objects in the sam¬ 
ple would have been detectable at A(M/)=2.1 mag, 2.5 
and 3.25 mag (85%, 90% and 95% of light absorbed lost). 
At 99% di mming, no obje c ts wo uld remain detectable. 

For the [Pierce fe Tullvi (119881) sample of disk galax¬ 
ies used bv lAnni si 1119991:. . one may derive a correspond¬ 
ing limit by considering the apparent magnitude limit at 
which i ncom plet eness e f fects s et in. Following the discus¬ 
sion in iPierce fe Tullvi (|1988T ), we adopt a completeness 
limit at to/ c >11, and argue that « 17 out of the 57 disk 
galaxies (i.e ~ 30%) used bv lAnnis! (119991) are sufficiently 
bright to have been included in the sample even if they 


had been dimmed by 1.5 mag. By using the photometric 
data fr om th e SFI++ compilation for the subset of the 
lAnnid (119991 ) disk galaxies (43 objects) that are included 
in SFI-I—(-, we can aslo perform a completion analysis 
based on surface brightness similar to that used above 
for our own sample. Th is sugg ests that a slightly larger 
fraction (« 40%) of the A lmisl (119991 ) disks may remain 
detectable after 1.5 mag of dimming. He nce, th e non - 
detection of Kill candidates reported by Anni 1 (fl999l) 
converts into a <4-6% upper limit on the prevalence of 
Kills in disk galaxies. 

This means that, due to the existence of the apparent 
outlier population shown in Fig. [1} our upper limit on the 
prevalence of star-fed Kill civilizations in disk galaxies is 
at most a factor of » 2 ti mes stronger than that inferred 
from the original [Annidl (1 9991 ) study (< 3% for our sam¬ 
ple compared to < 4-6% for the Annis sample), despite 
the fact that we are analyzing a disk galaxy sample that 
is larger by a factor of « 20. As discussed in Sect. 12.21 
many of the outliers we find are likely to be due to inac¬ 
curate distances, and this is an effe ct that will plague all 
large TF surveys - the 1 Annis (1999) approach, with focus 
on just a single, well-defined cluster, provides a cleaner 
sample in this respect. The presence of a small fraction of 
disk galaxies that may appear underluminous due to low 
star formation activity or unusually high dust content 
just compounds the problem. Applying the TF method 
to even larger samples may therefore not do much to 
improve the upper limits on the Kill civilizations, un¬ 
less auxiliary data allows for stronger constraints on the 
nature of the outliers. 

In Sect. [3] we discuss the use of morphological infor¬ 
mation and infrared data to further probe the Kill sta¬ 
tus of our TF outliers. As none of our 11 TF outliers 
stand out as strong Kill candidates based on these addi¬ 
tional tests, we may tentatively conclude that a fraction 
< 1/400 (i.e. < 0.3%) of disk galaxies host star-fed Kill 
civiliza tions ( i.e. a constraint 10-20 times stronger than 
that of lAnnislll999D . based on the assumption that none 
of the 11 candidates are bona fide Kill objects. 

3. THE NATURE OF THE KARDASHEV TYPE III GALAXY 
CANDIDATES 

As described in Sect. [2 a small fraction of local disk 
galaxies appear to be significantly underluminous (by 
a factor of four, i.e. by 1.5 mag) compared to the 
TF relation. Individual objects that are TF outliers 
in this res pect have been discussed in the literature 
before (e.g. iBarton et al.l 12001b IKannappan et al.l l2002t 
iPizagno et all I2007D . without much consensus on the 
likely cause of their low luminosities, although it has 
been noted that disks of morphological class Sa display 
a potential systematic o ffset in this direction of the TF 
(IKannappan et ahl 1 200 2D . As seen in Tabic Q] only 1— 
2 of the Kill candidates in our final sample bear this 
classification. In the following, we will discuss auxiliary 
observations to further probe the nature of our candi¬ 
dates. 

3.1. Morphological signatures 

As the colonization wave of a civilization in the transi¬ 
tion phase between KII and Kill sweeps across a galaxy, 
conspicuous morphological anomalies may potentially re¬ 
sult, but only as long as the Dyson spheres remain con- 
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TABLE 1 

Candidates for host galaxies of star-fed Kardashev type III civilizations 


No. 

Ra 

(J2000) 

Dec 

(J2000) 

Name 

log Wtf 
( km s -1 ) 

cr(log W T f) rn Ic 
l (km s -1 ) 

Mj 


^grp,M 

(km s -1 ' 

^^gal 

) (km s -1 

Type 

) 

Group 

Class 

1 

00 09 48.2 

+27 49 55 

NGC 0022 

2.696 

0.033 

12.42 

-20.64 

0.13 

2862 

7980 

Sb 

40702 

c 

2 

00 11 45.1 

+28 29 56 

UGC 00108 

2.652 

0.007 

13.16 

-19.90 

0.13 

2862 

7709 

Sb 

40702 

c 

3 

10 01 47.9 

+36 29 56 

UGC 05394 

2.137 

0.009 

14.45 

-17.14 

0.23 

1455 

1685 

Sc 

40032 

A 

4 

11 45 41.2 

-28 22 03 

ESO 440-G004 

2.360 

0.059 

13.11 

-18.98 

0.29 

1827 

2182 

Scd 

30590 

B 

5 

12 32 03.2 

+16 41 13 

NGC 4502 

2.258 

0.024 

12.89 

-18.17 

0.28 

1140 

1944 

Scd 

30654 

A 

6 

12 48 22.9 

+08 29 15 

NGC 4698 

2.755 

0.046 

9.20 

-21.86 

0.27 

1140 

1330 

Sa 

30654 

A 

7 

12 54 48.5 

+ 19 10 34 

IC 3877 

2.324 

0.012 

12.35 

-18.71 

0.28 

1140 

1216 

Sc 

30654 

A 

8 

13 26 12.8 

-27 29 06 

AGC 530433 

2.464 

0.023 

16.21 

-19.79 

0.09 

11080 

9992 

SO/a 

1736 

B 

9 

22 13 30.6 

-27 33 30 

ESO 467-G023 

2.655 

0.013 

11.82 

-20.05 

0.20 

1651 

4976 

Sb 

31043 

c 

10 

22 15 48.2 

-27 30 44 

ESO 467-G034 

2.380 

0.025 

13.87 

-18.00 

0.21 

1651 

4859 

Sb 

31043 

c 

11 

07 03 26.7 

-48 59 40 

AGC 470027 

2.560 

0.035 

15.32 

-20.91 

0.08 

12346 

12574 

Sb 

3407 

A 


The entries in the group column correspond to group identifiers from the SFIH—b catalog. 


fined to specific regions within that object. Once a 
galaxy is completely colonized, a uniformly dim disk is 
the more likely outcome. In principle, a completely dark 
disk may also arise (which would be undetectable with 
the current method), but this requires that every type of 
star is targeted for astroengineering, and that the Dyson 
spheres capture 100% of the optical light, which may be 
difficult to achieve in practice. 


In the context of Dyson spheres constructed by self- 
replicating von Neumann probes, colonization will typ¬ 
ically e ngulf a galaxy like the Milky Way in < 10 7 -10 8 


yr (e.g. Tinledl 198f* 

Valdes & Frcitaslll98Ct Barlow! 2013t 

INicholson & Forgan 

2013(1. resulting in a uniformly dim 


system on time scales much shorter than the current ages 
of large disks (~ 10 10 yr). Hence, we are unlikely to catch 
the system in the KII-KIII transition phase where mor¬ 
phological signatures would be present. 

For non-von Neumann type colonizers, estimates on 
the time required to colonize a galaxy vary greatly 
depending on assumptions such as vessel propulsion 
method, number of probes and pop ulatio n dynamics - 

from ~ 10 6 yr up to ~ _ IQ 1 0 vr (lHar tl 1 197,4 I.Tonesl 

19761 iNewman &; Saeanl fl981k iSaean fe Newmanl 19831 

Cotta fe Morales!l2009HForgan et alJl2013t iWright et al l 
2014al) . But even if the time scale for complete coloniza¬ 
tion is long, would this necessarily result in long-lived 
morphological anomalies during the colonization phase? 
In many cases, the answer is no. In the case of very 
slow-moving colonizing vessels, stellar dispersion would 
still spread the colonized st ars t hroughout most or all of 
the galaxy in ~ 10 8 -10 9 yr dWright et al.ll2014af) . A com¬ 
parison of this time scale with the ages of disks (~ 10 10 
yr) already indicates that < 10% of KII-KIII transition 
objects should display large-scale morpholgical anoma¬ 
lies. 

What if the coloinzation efforts were to stop abruptly, 
once a certain fraction of the stars had been converted 
into Dyson spheres? In the case of fast-moving vessels 
(moving faster than the stars within the galaxy), the 
colonization pattern at the stalling point could be highly 
irregular, but differential rotation and stellar dispersion 
would still carry the luminous stars into the colonized 
regions, and result in a uniformly dim disk in few times 
the rotational period (typically ~ 10 8 -10 9 yr). How¬ 
ever, there are special situations in which morphologi¬ 
cal curiosities could be long-lasting, due to the limited 
mixing between the dynamical subcomponents of galax¬ 


ies. As an extreme case, consider a colonization wave 
that spreads exclusively throughout the bulgt0 of a disk 
galaxy, but stops abruptly once some fraction of bulge 
stars have been turned into Dyson spheres. This would 
render the bulge significantly under luminous, and while 
some of these Dyson spheres may wander into the disk, 
they will not all end up on disk-like orbits for any ap¬ 
preciable amount of time. After all, the reason why sig¬ 
nificant age and metallcitiy differences are observed be¬ 
tween bulge and disk po pulations in individual galaxies 
(e.g iSancliez-Blazauez et ahll2014f) is that stars in these 
subcomponents have retained their orbital identities for 
several billions of years. 

In Fig. [2j we use a simple numerical model of galaxy 
colonization to illustrate a long-lasting morphological 
anomaly that arises because of a stalling colonization 
wave. In these simulations, the stellar population of a 
galaxy is represented by 5 x 10 5 star particles, each rep¬ 
resenting a large number of stars (~ 10 6 ). The star parti¬ 
cles are distributed throughout the thin disk, thick disk, 
bulge and halo of a galaxy with dimensions similar to 
the Milky Way (disk scale lengths 3 kpc; thin disk scale 
height 0.3 kpc; thick disk scale height 1 kpc; bulge radius 
1.5 kpc, halo radius 30 kpc), and assigned velocities de¬ 
pending on which component they belong to. The masses 
of the bulge, thick disk and stellar halo are assumed to 
be 30%, 5% and 1% of the mass of the thin disk, and 
the number of particles in these components scale ac¬ 
cordingly. The particles within the disks are assigned ro¬ 
tational velocities giv en by the universal rota tion curve 
for spirals specified in lPersic fe Saluccl (119951) . An addi¬ 
tional velocity component in the polar direction is also 
assigned to make each particle oscillates around its ini¬ 
tial position with a maximum polar velocity of 30 km 
s^ 1 . A random radial velocity component with velocity 
dispersion 10 km s^ 1 is also included in every time step. 
The bulge and halo particles are assigned two velocity 
components: a fixed rotational velocity of 100 km s^ 1 
and 200 km s^ 1 respectively in randomly-inclined orbits, 
plus the same random radial migration velocities (with 
respect to their inclined orbital planes) as disk particles. 

A random star particle at 8 kpc from the centre (the 
distance of the Sun from the centre of the Milky Way) 
is chosen as the starting point of the colonization, and 

5 Here taken to mean a classical bulge formed through mergers 
or monolithical collapse, as opposed to a pseudo-bulge which forms 
from secular processes in the disk. 
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the evolution of the colonization wave is then tracked 
throughout subsequent timesteps as probes and star par¬ 
ticles move about. We adopt a colonization strategy in 
which each new colony immediately sends out a very 
large number of probes with velocities 10~ 3 c (i.e. some¬ 
what faster than the average star particle) in all direction 
as soon as it is captured. This scheme allows 50% of the 
star particles in the model galaxy to turn into colonies 
in less than 10 8 yr, after which we assume further colo¬ 
nization attempts to cease. This 50% level is for plotting 
purposes chosen as the end stage of the colonization ef¬ 
forts, but the 75% case that more closely matches our 
observational search criteria would not produce results 
that are fundamentally different. 

The model galaxy is shown in face-on and edge-on ori¬ 
entations throughout Fig. [2ji-f, in various stages before 
and after the stalling phase. The colour scheme for the 
star particles indicates the projected surface mass density 
at their radius from the centre, with purple/blue indicat¬ 
ing the highest stellar density, and red the lowest. Once 
a star particle is tagged as colonized, it is no longer con¬ 
sidered to contribute to the integrated optical luminosity 
of the galaxy, and the colour of that star particle is then 
set to black. Fig. Ek depicts the disk just 5 x 10 6 yr (at 
< 1% colonization) after the release of the first coloniz¬ 
ing vessels. Fig.Eb shows the corresponding distribution 
of colonies at 2 x 10 7 yr (s=s 15% colonization). Since the 
adopted scheme allows for more efficient colonization in 
regions of high stellar density, the colonization wave is 
making its way towards the inner regions of the galaxy 
and has at this point begun overtaking the bulge. The 
inner regions of the galaxy are then relatively quickly 
turned into stellar real estate, making the galaxy reach 
50% colonization just a few Myr later (Fig. Et)- We 
assume that no more star particles are converted into 
Dyson spheres beyond that point, and show the state of 
the same simulated galaxy 2.5 x 10 8 yr, 5 x 10 s yr and 
1 x 10 9 yr later throughout Fig. EJi-f. Differential rota¬ 
tion pulls the colonized patch of the disks into a dark 
spiral, which - when coupled to random radial migration 
of stars within the disk- will eventually turn into a uni¬ 
formly faint disk. The situation at 2.5 x 10 8 yr after the 
stalling point (Fig. Eli) already closely resembles a system 
of galactic dust lanes (although without the associated 
optical reddening), and would not necessarily stand out 
as anomalous upon morphological inspection of an image 
in a single filter. Since the simplified simulations pre¬ 
sented here do not take density waves, bar instabilities 
or perturbations due to satellite galaxies into account, it 
could be argued that the “dark spiral” of colonized stars 
may disperse on a time scale somewhat shorter than sug¬ 
gested by Fig. El However, if a new bulge population 
does not form (due to mergers or secular evolution of the 
disk), the central region will remain under luminous for 
billions of years, since the mixing of particles in the disk 
and bulge is empirically known to be a very slow process. 

To summarize, morphological curiosities are possible, 
albeit not necessary, signatures for Kill galaxies. Based 
on optical images of our Kill candidates publicly avail¬ 
able through the NASA/IPAC Extragalactic Database 
(NED), none of them appear to display any conspic¬ 
uous morphological anomalies. Object 4 (ESO 440- 
G004) in Table |T| has a disturbed/irregular appearance, 
with sprawling spiral arms connected to a central bar 


ijMatthews fe Gallagherlll997ll . but this is more likely to 
stem from a merger and does not resemble anything in 
our colonization simulations. Hence, we find no reason 
to attribute the optical morphologies of any of our Kill 
candidates to galaxy-scale astroengineering projects. 

3.2. Infrared excess 

While Dyson spheres may be very efficient at absorb¬ 
ing the ultraviolet and optical radiation of stars, waste 
heat would probably still need to be radiated away. For 
Dyson spheres constructed out of currently known ma¬ 
terials, the structure would need to operate at a tem¬ 
perature that places th e peak of the waste radiation at 
infrared wavelenghts. iDvsonl (|1960H suggested that one 
should search for excess radiation at ps 10 /im, which 
corresponds to the radiation peak of T « 300 K black 
body, and a number of searches have already used the 
InfraRed Astronomical Satellite (IRAS, operating at 12- 
100 fi m) to search for Dyson spheres operating at ~ 50- 
600 K (iSlvslill 19851 lTimofeevll2000t I Jueaku fe Nishimural 
12004 iCarriaaiJ 120091) . The G-H AT survey (Glimpsi ng 
Heat from Alien T echonoloqies\ iWright et al.1 l2014allb[ : 
iGriffith et alJl2015l ) is currently also looking for such in¬ 
frared excess signatures (from both KII and Kill civi- 
liziations) using data from the Wide-held Infrared Ex¬ 
plorer (WISE). 

In the case of galaxies, predicting the exact amount 
of infrared excess at a specific wavelength turns into a 
highly non-trivial spectral synthesis problem, since the 
flux excess depends both on the type of stars targeted by 
Dysonian astroengineering, on the temperature at which 
the Dyson sphere operates and - to some extent - on the 
intrinsic dust properties of the galaxy. Our Kill can¬ 
didates are underluminous by > 1.5 magnitudes in the 
/-band, but unless Dyson spheres are constructed around 
all types of stars, this provides but a poor handle on the 
bolometric luminosity absorbed by the Dyson spheres 
(and re-emitted in the infrared). As an example, con¬ 
sider the case where the type of stars that contribute the 
most to the /-band flux are exclusively targeted by Dyson 
sphere projects. This would correspond to an absorbed 
total luminosity lower than that inferred from reducing 
the integrated, bolometric stellar luminosity of the entire 
galaxy by 1.5 mag from its original value. If, on the other 
hand, the stars dominating the /-band are among the 
least preferred targets of astroengineering, the absorbed 
fraction may correspond to more than a 1.5 mag reduc¬ 
tion of the bolometric, integrated stellar luminosity of 
the galaxy. Furthermore, part of the optical flux of star¬ 
forming galaxies comes in the form of nebul ar emissio n 
rathe r than direct starlight (e.g. iZackrisson et al.ll200ll . 
2008). If the hot stars (spectral types O and B) that con¬ 
tribute the most to the photoionization of the interstellar 
medium are enshrouded in Dyson spheres, the total flux 
decrease at optical wavelengths will be greater than that 
inferred by summing up the direct optical flux contribu¬ 
tion from these stars. Since much of the dust heating 
responsible for the mundane mid/far-infrared fluxes of 
disk galaxies would also decrease as these stars are con¬ 
verted into Dyson spheres, the dust contribution to the 
mid/far-IR fluxes of Kill galaxies would also be lowered. 

Here, we will not attempt a full-blown spectral syn¬ 
thesis analysis of our Kill candidates, but - given that 
a bona fide Kill candidate would be expected to dis- 
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Fig. 2.— Snapshots from a simulation in which star particles within a disk galaxy are colonized by rapidly replicating probes moving 
at 10 -3 c, but where colonization efforts cease once 50% of the star particles have been colonized. The purple-red color scale indicates 
the stellar density in the uncolonized disk, with purple marking the highest density and red the lowest. Colonized particles are marked 
in black, a) The disk in face-on and edge-on orientations 5 Myr after the release of the colonization wave from a position 8 kpc from the 
centre, b) The state of the galaxy at ^ 15% colonization, after 20 Myr. c) The state of the galaxy at pa 50% colonization, after 25 Myr, 
at which point no more star particles are assumed to be colonized. In subsequent panels, the state of the simulated galaxy is shown at 
2.5 x 10 8 yr (d), 5 x 10 8 yr (e) and 1 x 10 9 yr (f) after this stalling point. The time stamp in each panel indicates the time since the launch 
of the first colonizing vessels. While differential rotation and radial migration make the colonized stars in the disk tend towards a uniform 
disk distribution throughout panels (e)—(f), the bulge colonies remain concentrated to the centre, leaving the inner galaxy significantly 
underluminous at optical wavelengths for billions of years. See main text for further details. 
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Fig. 2.— Continued 
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TABLE 2 

IR-TO-OPTICAL LUMINOSITY RATIOS OF THE Kill CANDIDATE 
GALAXIES FROM TABLE [T| 


No. 

L^.q/Li 

L 22 /L 1 

1 

0.94 

67 

2 

1.02 

43 

3 

0.24 

62 

4 

0.22 

46 

5 

0.35 

48 

6 

0.43 

2.7 

7 

0.11 

- 

8 

2.7 

- 

9 

1.9 

690 

10 

0.54 

91 

11 

4.9 

900 


Missing entries in the L 22 /Lj column indicate non-detections 
( S/N < 2) in the 22 (im band. 

play some sort of IR excess - simply investigate whether 
there is something clearly anomalous about the infrared 
properties of our Kill candidates. While IRAS, Spitzer 
and ISO data are available for some of our candidates, 
all of them have been observed by WISE at 3.4, 4.6, 
12 and 22 pm (although a few are undetected in the 12 
and 22 pm passbands). Here we will focus on the 4.5 
and 22 pm bands, since these are the least affected by 
emission features from polycycli c aro matic hydrocarbon 
(PAH) molecules (iWright et alJ2010h . and should be sen¬ 
sitive to infrared excesses associated with Dyson spheres 
operating at T « 80-1200 K under the assumption that 
these may be treated as single-temperature black bod¬ 
ies. We note, however, that this assumption may break 
down in cases where there are variations in the charac¬ 
teristic Dyson sphere temperatures among the colonized 
stars within a galaxy. Single-temperature black bodies 
may also be poor representations of stars enshrouded by 
multiple Dyson shells with partial coverage, where each 
layer is operating at a different temperaturqj, or so-called 
stellar engines with temperature differences bet ween the 
inner and outer surfaces of partial Dys on shells (|Badescul 
ITM! IBadescu fo Cathcartll2000L l2006h . 

In Table[2] we present the L^/Li and L 22 /L 1 in-band 
luminosity ratios (in solar units) for our Kill candidates: 

£ 4 . 6/22 lO _O - 4 ( m 4.ii/22-Af4. 6/22i0 ) 

Li ~ IQ 

w here m j c , is the /-band apparent magnitudes from 
iSpringob et al.l (|2007i ) corrected for Galactic extinction 
and inclination (internal dust attenuation), 7774.6 and in 22 
represent the 4.6 and 22 pm apparent magnitudes from 
WISE whereas A// and A/j , 0 are the corresponding /- 
band absolute magnitudes of the target galaxy and of 
the Sun in these passbands (A//, 0 = 4.10, A/4.6,0 = 3.27 
and A/22,0 = 3.25 in the Vega system). Missing entries 
in the £ 22 /-// column for two of the Kill candidates in¬ 
dicate that these are undetected (S/N < 2) at 22 pm. 

As Dyson spheres are expected to lower the I -band lu¬ 
minosities and boost the mid/far-infrared luminosities of 
stars, one naively expects the L^^/Li and/or the £22/ £/ 
ratio of a galaxy to grow substantially when a large frac¬ 
tion of its stars are subject this form of astroengineering. 
Since these ratios are independent of distance, they are 

6 As in the case of Robert Bradbury’s “Matrioshka brain” 


not affected by any lingering distance errors in the TF 
analysis (Sect. 12.2|) . One caveat is that if the Dyson 
spheres operate at a very low temperature (< 70 K), the 
excess heat radiation would mostly affect wavelengths 
longward of 22 pm, and whatever flux boost these low- 
temperature Dyson spheres may contribute at 22 pm 
could in principle be offset by a drop in the mundane 
22 pm dust radiation as stars that under normal circum¬ 
stances keep this dust heated are enshrouded in Dyson 
spheres. The second problem is that the intrinsic scatter 
among disks in £m /£ n nt.i rai luminosity ratios a mounts 
to a factor of ~ 100 fe.g. iKennicutt et all 120031 ) . which 
means that - without detailed spectral synthesis mod¬ 
elling - a truly extreme IR excess is required to stand 
out. 

As seen in Tabled Li^/Li varies by a factor of « 40 
and L 22 /L 1 by a factor of « 300 among our candi¬ 
dates, with the most extreme objects reaching L^^/Li Ri 
5 £ 4 . 6,0 /£j, 0 and £ 22/£ 7 ~ 9 00 £ 22 . 61 /' £r. 0 - However, a 
similar analysis of lSpringob et all ( 2007fl disks with more 
typical /-band luminosities with respect to the TF rela¬ 
tion indicates that such objects occasionally also attain 
L^.q/Li and L 22 IL 1 ratios this high. 

The Kolmogorov-Smirnoff test rejects (at 5% signif¬ 
icance) the hypothesis that the L\^/Li and L 22 /L 1 
distributions of our Kill candidates is drawn from the 
same par ent population as a random sample of disks 
from the ISpringob et all (120071 1 catalog. However, this 
is primarily due to an excess of objects with very low 
IR-to-optical flux ratios in our Kill candidate sample 
(Li.e/Li < 1 £ 4 . 6 , 0 /£ 7 , 0 , L 22 /L 1 < 100 £ 22 , 0 /£/,©)> 
i.e. objects with ratios of a kind not expected from bona 
fide Kill host galaxies. 

To explore whether the mid-IR spectra of our Kill can¬ 
didates shows some other interesting anomalies, we have 
plotted these objects into the TO 3.4 — 7714.6 vs. 7774.6 — 77742 
color-color di agram ( Fig. EH), with a few of the diagnostic 
regions from IWright et all (120101) superposed. Due to a 
lack of reliable 12 pm data for one (object 3), only 10 
objects are included in this figure. 

Most of these Kill candidates fall in the region oc¬ 
cupied by disk galaxies, which indicates that their mid- 
IR flux ratios in these WISE passbands are normal for 
objects of this type. However, one of the Kill candi¬ 
dates (object 6) is located close to the region typically 
inhabited by elliptical galaxies, where stellar continuum 
rather than dust emission dominate the WISE fluxes. 
Since this object has a fairly low L\^/Li ratio and a 
very low £ 22 /£/ ratio, it is not a strong Kill candidate 
based on its IR properties. Instead, a possible explana¬ 
tion for both its anomalously low /-band luminosity and 
lack of pronounced dust emission is that star formation 
has been quenched in this disk galaxy, thereby effectively 
shutting off dust heating by massive stars and allowing 
the Rayleigh-Jeans tail of the stellar radiation to domi¬ 
nate the WISE bands. In deed, t his i nterpretatio n is sup - 
ported by the analysis by iKoopmann fe Kennevl (1200411 , 
which classify it as “anemic” based on the extremely low 
equivalent width of its Ha emission line. Three of the 
candidates (objects 9, 10 and 11) are located in the part 
of the spiral galaxy region also occupied by Luminous In¬ 
frared Galaxies (LIRGs), and two of these (objects 9 and 
11) are also the ones with the highest £ 22 /£/ ratios. Out 
of these, only object 9 has IR data at longer wavelengths 
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(IRAS; 12-100 /im). This candidate is not a LIRG ac¬ 
cording to the standard defintion (L IR > 10 11 L 0 ), but 
it does have a very high IR-to-optical luminosity ratio for 
its type (total IR to B-band ratio Ttr/Db ~ 5 com pared 
to ss 0.4 for a typical Sb galaxy; Ide Jong et al.1 [19841. 
However, o ther disk galaxies with equally extreme ratios 
are known (iKennicutt, et, al.lfeOO.'ili . 

Despite the complications in predicting the exact SEDs 
of Kill galaxies, one may still set lower limits on the 
expected L^q/Li and T 22 /E/ ratios by adopting the very 
conservative assumption that only the /-band luminosity 
is absorbed by Dyson spheres, and that this is re-radiated 
as a single-temperature black body. In the case of black 
body Dyson spheres of temperature T = 75, 100, 300 
600 K, a 1.5 mag decrement in the /-band results in 
minimum luminosity ratios of T 22 /A 7 > 300, 800, 850 
and 210. For T = 400, 600 and 1000 K black bodies, 
L/i.§/Li > 4.7, 120 and 140, respectively. These lower 
limits already rule out most of the objects in Table [2] as 
Kill candidates with single-temperature Dyson spheres 
in the T > 75 K range. A detailed comparison reveals 
that T w 400-1200 K Dyson spheres are ruled out for all 
candidates based on the L^, e /Lj ratios, and that 75-800 
K Dyson spheres are ruled out for all but objects 9 and 
11 based on the L 22 /L 1 ratios. 

The data in Table [2] would make object 9 consistent 
with Kill activity under the assumption of Dyson spheres 
operating at very low temperatures (T < 100 K). How¬ 
ever, since this object belongs to class ’C’ due to a du¬ 
bious distance estimate (see Sect. H, we do not consid¬ 
ered it a strong Kill candidate. This leaves object 11, 
which has the highest L 22 /Lj and L^/Lj ratios among 
our Kill candidates (consistent with Kill activity in the 
case of T < 400 K Dyson sphere, under the very conser¬ 
vative assumption that only the I-band flux is absorbed) 
and has a redshift without any siginficant velocity off¬ 
set from that of its the parent group. However, the fact 
that this object straddles the Kill candidate delimiter 
line (Fig. |I]) makes its candidate status sensitive to the 
details of how the extinction corrections are made. Our 
luminosity-dependent extinction corrections differ some¬ 
what from those tabulated in the SFI++ catalog (due to 
slight differences in the way distances are assigned), and 
while this difference is in general very small (on average 
~ 0.02 mag throughout our sample and insignificant for 
the Kill candidate status of the other 10 objects in Ta¬ 
ble [lj. the discrepancy for this particular object happens 
to be to sufficient to push it ~ 0.3 mag closer to the TF 
if the original SFI++ correction is used. While hardly 
a strong Kill candidate, this is still the most interesting 
object uncovered in our search, and could be an interest¬ 
ing target for follow-up observations. 

4. DISCUSSION 

Without knowing what an extremely advanced civi¬ 
lization would use the energy collected from large num¬ 
bers of Dyson spheres for, it is difficult to judge whether 
this type of astroengineering project would be allowed 
to grow indefinitely, or cease once some desired energy 
requirements have been met. If the pu rpose is t o har¬ 
vest energy to drive supercomputers ([Sandberg 119991 
lBradburvll2001f) the long communication times between 
distant Dyson spheres may eventually make further ex¬ 
pansion of this type of Dysonian colonization unattrac- 



Fig. 3. — WISE m 3 . 4 — m. 4.6 vs. 7714.6 — mi 2 color-color dia¬ 
gram for the 10 of our Kill candidates with WISE detections in 
these passbands. The objects have been labeled according to their 
numbers in Table [I] The colored patches mark the regions typi¬ 
cally inhabited by elli pticals, disks a nd Luminous Infrared Galaxies 
(LIRGs) according to lWright et al.l <120100 . In some cases, the error 
bars on the measured colors are smaller than the symbols used. 

tive. 

Civilizations that cease their colonization efforts early 
in the transition stage from Kar dashev type II to III 
would not be singled out using the lAnni 1 (fl999 1 Kill cri¬ 
terion (which requires that 75% of the optical light has 
been absorbed), but could potentially_still be detectable 
in nearby, well-resolved galaxies. iCarriganl (l2010i l ar¬ 
gue that regions within such galaxies where large num¬ 
bers of stars have been converted into Dyson spheres 
would appear as dark patches in hig h-resolution opti¬ 
cal images. As discussion in Sect. 13.11 large (kpc-scale) 
patches of this type would only survive over cosmological 
time scales under very special circumstances. However, 
a gravitationally bound star cluster subject to astroengi¬ 
neering could retain its Dysonian signatures for billions of 
years. Even though optical dimming coupled to an IR ex¬ 
cess superficially mimics the effects of dust attenuation, 
Dysonian astroengineering projects may still give rise 
to quantitative anomalies in integrated UV/optical/IR 
colors or spectra of such objects, that would be differ¬ 
ent from those produced by standard dust attenuation 
curves. This is especially true if only a subset of stellar 
types are targeted by astroengineering projects. 

If Dysonian astroengineering ambitions remain firm 
over time and are not hindered by scale, the colonization 
wave may in princip le spread beyo nd t he b ou ndar ies of 
its galaxy of origin. lArmstrong fe Sandberd (1201311 out¬ 
line various schemes through which von Neumann probes 
may spread to large number of galaxies over cosmolog¬ 
ical time scales. This provides another opportunity to 
confirm the status of Kill candidates - if neighbour¬ 
ing galaxies would also appear anomalously dim, this 
would constitute supporting evidence for astroengineer¬ 
ing. Promising methods to investigate this may be to 
study the luminosity functions of dwarf galaxies in the 
vicinity of Kill candidates. In cases where the Kill can¬ 
didate happens to be located close to another large disk 
or elliptical gal axy, a Tully-Fisher or fundamental plane 
analysis (|Annislll999tl would be very interesting. Indeed, 
several of the Kill candidates in Table [T] belong to the 
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same galaxy group. Among objects in the ’A’ category 
(indicating that the group has > 4 members in the cat¬ 
alog and therefore a more reliable TF distance to the 
group centre), this is the case for objects 5-7. However, 
these objects also happen to belong to the most nearby 
group featured in Tabic [Tj (among th e 3% o f groups with 
the lowest inferred distances in the lSpringob et all 120071 
catalog), which makes the distance of individual galaxies 
from the group centre a significant error source. With¬ 
out independent (i.e. non TF-based) distance estimates 
to the individual galaxies in the groups containing sev¬ 
eral Kill candidates, we are unable to assess how distant 
these galaxies really are from one another. 

5. SUMMARY 

By searching a Tully-Fisher sample of 1359 disk galax¬ 
ies for the Dysonian astroengineering signatures expected 
from Kardasliev type III civilizations, we are able to set 
a conservative upper lim it of < 3 % o n the fraction of 
local disks that meet the lAnnisl (11991 criteria for such 
civilizations. In this sample, a total of 11 objects are 
found to be significantly underluminous (by a factor of 
4 in the /-band) compared to the Tully-Fisher relation, 
and therefore qualify as Kardasliev type III host galaxy 
candidates according to this test. However, by scruti¬ 
nizing the optical morphologies and WISE 3.4-22 /im 
infrared fluxes of these objects, we find nothing that 


strongly supports the astroengineering interpretation of 
their unusually low optical luminosities. Hence, we con¬ 
clude that their apparent positions in the Tully-Fisher di¬ 
agram likely have mundane causes, with underestimated 
distances being the most probable explanation for most 
of the candidates. Under the assumption that none of 
them are bona fide Kill objects, we set a tentative upper 
limit of < 0.3% on the fraction of disk galaxies harbour¬ 
ing Kill civilizations. 
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